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ABSTRACT 

Context. The study of variations in total solar irradiance (TSI) is important for understanding how the Sun affects the Earth's climate. 
Aims. Full-disk continuum images and magnetograms are now available for three full solar cycles. We investigate how modelled 
TSI compares with direct observations by building a consistent modelled TSI dataset. The model, based only on changes in the 
photospheric magnetic flux can then be tested on rotational, cyclical and secular timescales. 

Methods. We use Kitt Peak and SoHO/MDI continuum images and magnetograms in the SATIRE-S model to reconstruct TSI over 
cycles 21-23. To maximise independence from TSI composites, SORCE/TIM TSI data are used to fix the one free parameter of the 
model. We compare and combine the separate data sources for the model to estimate an uncertainty on the reconstruction and prevent 
any additional free parameters entering the model. 

Results. The reconstruction supports the PMOD composite as being the best historical record of TSI observations, although on 
timescales of the solar rotation the IRMB composite provides somewhat better agreement. Further to this, the model is able to account 
for 92% of TSI variations from 1978 to 2009 in the PMOD composite and over 96% during cycle 23. The reconstruction also displays 
an inter-cycle, secular decline of 0.20*?'™ WnT 2 between cycle 23 minima, in agreement with the PMOD composite. 
Conclusions. SATIRE-S is able to recreate TSI observations on all timescales of a day and longer over 31 years from 1978. This is 
strong evidence that changes in photospheric magnetic flux alone are responsible for almost all solar irradiance variations over the 
last three solar cycles. 
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1. Introduction 

The Sun is the primary source of radiation entering the Earth's 
climate system, reaching to all levels of the atmosphere and into 
the ocean ( |Haigh|1994HBond et al.|200T]|Labitzke et al.|2002) . 



While much of the modulation from the Sun seen in the Earth's 
atmosphere is a result of variation in spectral solar irradiance 
(SSI), total solar irradiance (TSI) encompasses the combined 
variation across all wavelengths and thus is a measure of the 
total forcing. Consequently, TSI is an important quantity in cli- 
mate studies. 

TSI is also easier to measure than SSI and should, in prin- 
ciple, provide a more stable and accurate record of solar vari- 
ability. However, while direct, above-atmosphere, radiometric 
observations exist from 1978, continuous observations have re- 
quired the use of several satellites. This has lead to an inevitable 
need to account for different absolute level calibrations, degra- 
dation corrections and switch-off periods in each instrument. 
Accounting for each of these problems naturally leads to some 
interpretation and, as a result, three competing TSI compos- 



ites have been proposed (Frohlich & Lean||1998 


|Mordvinov|2003 


Dewitteetal.|2004i. 


Wenzler et al. 


(20091 and Krivova et al.|(|2009) 



the best record of TSI, but this required fitting a model to each 
composite separately and was not independent of the compos- 



ites. Also, the reconstructions used in these investigations ended 
in 2003 and hence do not cover the unusual minimum between 
cycles 23 and 24. 

Many studies have shown that the Sun's surface magnetic 
flux, in the form of sunspots and faculae, can explain TSI varia- 
tion on rotational time scales of days and weeks ( Willson et al. 
[T98"T||Foukal &Lean| 19861 |Fligge et al.|2000|> and cyclical pe- 
riods covering the solar cycle (Lean et al.|| 1 998 ; Wenzler et al. 



2006 ). All three TSI composites suggest that the Sun undergoes 
some measure of inter-cycle, secular variability, but all three are 
in disagreement as to the magnitude and direction of this change. 
Some indices, e.g. open solar flux ( [Lockwood et al.|2010 1, show 
that the minimum of 2008 saw the largest significant secular 
change of the satellite era and this is reflected in the PMOD 
composite by a decline in TSI relative to the minima of 1996 
and 1986. Although the decline in T SI between 199 6 and 2008 
is small, estimated to be ~0.2 Wm~ 2 ( |Frohhch|2009| , the impor- 
tance of this small change lies in the accumulated change over 
centuries. It has been suggested that the Sun may be entering a 
grand minimum state, though there is still some debate on this 
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( Feynman & Ruzmaikin|20TT] Solanki & Krivova|201 1 1. 

To understand the impact of the Sun on climate change, it is 
imperative that the physical mechanism behind the variation of 
the Sun on secular timescales is understood to a sufficient level 
to be able to confidently reconstruct solar variability into the 
past. This would then allow reconstructions that extrapolate TSI 
proxy data, such as sunspot number and the geomagnetic aa in- 
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dex, to be refined and improved ([Ta pping et al. 2007 ; Steinhilber 
eTatl|2"0"09{ |Vieira & Solanki||2010| ). The potential range for a 
secular change in TSI over the last four hundred years, since the 
Maunder Minimum, could be significant in climate forcing and 
a better understanding of the Sun during the satellite era could 
lead to more precise estimates of the Sun's true effect on the 
Earth over that period, although there are diverging views on 
this. On the one hand, Shapiro et al. ( 201 1[ ) argue that because 
the Sun has been in a particularly active phase during the last 
century, one cannot currently extrapolate the full range of solar 



variability on secular scales. On the other hand, Schrijver et al 



( |201 1 ) suggest that solar activity during the 2008 minimum has 
reached levels similar to the Maunder Minimum, a period during 
which sunspots were conspicuous by their absence, and cannot 
go much lower. 

Explanations for secular trends other than solar surface mag- 
netic activity have also been proposed. E.g., Rozelot et al. ( 2004 ) 
suggested that changes in the Sun's radius could account for 
long-term changes in the Sun. Harder et al. (2009 1 proposed that 
there may have been a change in the photospheric temperature 
gradient that produced a spectrally dependent change, increas- 
ing irradiance at some wavelengths and decreasing at others. 
A global temperature change has also been proposed (Tapping 



et al. 2007 ) with Frohlich ( 2009 1 concluding that the recent min- 
imum was caused by a global temperature decline of 0.2 K in the 
Sun's effective temperature. 

The SATIRE-S model, denoted with an 'S' for the satellite 



era (see |Krivova et al. ( 201 lb\ to distinguish between the vari- 
ous SATIRE types), is a well established model that has proven 
very successful. Using Kitt Peak full-disk data, between 1974 
and 2003, the model has been able to reproduce 83% of TSI vari- 



2006 ). Using SoHO data, Krivova et al. 



2011) were able to account for 92% and 



ability (Wenzler et al. 
(2003) and |BaU et alj 
97% during the ascending and descending phases of cycle 23, 
respectively. 

|Steinhilber (2010) used recently recalibrated MDI synoptic 
data within a model framework similar to that of the SATIRE 
model over the full period of cycle 23, but was unable to repro- 
duce the recent decline seen by the PMOD composite. The cutoff 
threshold for low magnetic flux used was higher than previous 
reconstructions and could have removed a significant portion of 
the relatively weak magnetic flux in the network that might be 
responsible for secular variation. 

In this study, we combine all available full disk continuum 
images and magnetograms from Kitt Peak (KP) and the Solar 
and Heliospheric Observatory (SoHO) to compute TSI for a 
period of 36 years, extending the reconstruction, from 1974 to 
2009. 

In section|2]we briefly outline the SATIRE-S model. In sec- 
tion [3] we describe the full-disk images used in the model, the 
TSI data used to set up the reconstruction and the TSI compos- 
ite data used to make comparisons with. Continuum images and 
magnetograms from different sources each have their own intri- 
cacies that need to be accounted for before they can be com- 
bined. We explain how this is done in section|4] The datasets are 
combined to produce a consistent three-cycle reconstruction in 
section[5] Comparisons are made between the reconstruction and 
TSI composites in section|6]and discussion of the results and the 
conclusions are presented in section [7] 



2. Modelling Solar Irradiance using SATIRE-S 

Solar irradiance is modelled using the SATIRE-S model that as- 
sumes that only changes in surface magnetic flux are responsi- 



ble for changes in solar irradiance on time scales greater than 
one day.|Fligge et al.|(|2TXXfy|Krivova et al.| ( p003) and | Wenzler 
leTaL] d2004| |2005| |2006| ) explain the model in depth and |Ball 
|et al.| ( |2011| ) cover recent modifications. A detailed description 
can be found in these papers and references therein; here we give 
a short description only. 

SATIRE-S is a four-component model that uses daily con- 
tinuum intensity images and magnetograms. After correcting for 
limb darkening, sunspot penumbra and umbra are identified by 
contrast level in continuum intensity images. Faculae are identi- 
fied in magnetograms, where no spot pixels are located, above a 
threshold determined by the noise level, = 3cr noise (see sec- 
tions 3.1 and 4.5 1. Note that a line-of-sight correction is applied 
to the magnetograms, assuming that observed flux is radial in 
nature. This can produce artifacts at the limb, so all pixels with a 
limb angle of /i = cos 9 < 0.1, where 6 is the heliocentric angle, 
are ignored. All remaining pixels are considered to be the quiet 
Sun. 

Faculae are composed of magnetic elements too small to be 
resolved in the full disk magnetograms, so it is assumed that 
the facular filling factor increases linearly with increasing mag- 
netic flux up to a saturation point, B sat , above which the pixel 
is assumed to be 'full' and the filling factor is held at unity (cf. 
Fligge et al.| 2000). B sat is the only freely varying parameter 



in the model and is used to achieve good agreement with TSI 
observations. Magnetogram pixels of high magnetic flux are in- 
correctly assigned as faculae, so a cutoff, B cut , was introduced 
( 201 1[ ). This fixed parameter is discussed in see- 



in IBall et al. 



tionl431 

Intensities for each component, depending on // and on the 
wavelength A, are determined and a model atmosphere is em- 
ployed for each respective component as set out in |Unruh et al.| 
( 1999). The resulting spectrum is integrated on each day to pro- 



duce relative variations in TSI. 

3. Data 

3.1. Data used in SATIRE-S 

We use full-disk continuum intensity images and magnetograms 
employed in SATIRE-S, taken from two locations: the National 
Solar Observatory Kitt Peak Vacuum Tower (KP) and the Solar 
and Heliospheric Observatory (SoHO). KP data were taken us- 
ing the 512-channel diode array (KP/512) from 1974. KP/512 
was replaced by the spectromagnetograph (KP/SPM) which op- 
erated between 1993 and 2003. The Michelson Doppler Imager 
(MDI) operated onboard the SoHO satellite between 1996 and 
2011. 

3.1.1. SoHO/MDI 



MDI ( |Scherrer et IT] [1995) records 1024 x 1024 pixel full- 
disk filtergrams and magnetograms in the Nil absorption line 
at 676.78 nm with 4" resolution. Simultaneous continuum inten- 
sity image and magnetogram data are usually not available to 
download. Therefore, continuum intensity images taken within 
12 hours of the magnetograms are chosen to avoid the introduc- 
tion of large differences resulting from surface evolution or from 
rotation off the visible disk. The continuum intensity images are 
then rotated to align features in both images. Continuum inten- 
sity images have a measurement uncertainty of 0.3% and only 
a single image is needed to identify the sunspot components. 
For the magnetograms, we attain a reduced noise level by us- 
ing 5 -minute averaged magnetograms. The calibrated 1 -minute 
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magnetograms have a reported uncertainty of ~30 G; averaging 
over five-minutes reduces the uncertainty to ~15 G (these are 
recalibrated values following (Tran et al. 2005); see below). 

In this study, images are used from almost the entire work- 
ing period of the MDI instrument from 1996 May 20 to 2009 
October 31, using 4140 image pairs. 

Previous studies using SATIRE-S have employed level 1.5 
magnetograms (Fli gge et al.|2000||Krivova et al.|2003 Wenzler 
et al.|2004) . Recalibrations of the data have now produced the 
level 1.8.2 magnetogram data used in |Ball et"aL ( 201 l| l, Krivova 
et al.|(|201 la i and this work. The most significant change came 
after |Tran et al.| ( |2"005| ) compared the sensitivity of MDI magne- 
tograms with those of the Mount Wilson Observatory. This work 
led to a position-dependent correctior^jthat has increased the es- 
timated observed magnetic flux by a factor of ~1.6, necessitating 
a change in SATIRE-S parameters, notably B sat . 

SoHO was subject to a number of operational problems that 
MDI was not exempt from. It appears that these, along with some 
operational changes to the instrument, have resulted in changes 
to the magnetograms over time. We discuss the causes and ef- 



fects of these changes in section 4.5 



3.1.2. KP/SPM 



KP/SPM ( jJones et al.|[T992| operated between 1992 and 2003 
taking 1788 X 1788 pixel images in the Fe I line at 868.8 nm 
with a 1.14" pixel size. Continuum intensity images and magne- 
tograms were extracted from slit spectra that scanned across the 
solar disk in one hour. |Jones et"aT] ([T992 ) and We nzler et al.| 
(2004) both determined a background noise level of ~5 G ^ 
o"noise,sPM- Note that terms with the subscripts MDI, SPM and 
512 mean they are set for reconstructions using only that spe- 
cific dataset. 

Wenzler et al. ( |2005 i identified 2055 sets of images between 



1992 November 21 and 2003 September 21 that were relatively 
free of atmospheric effects. We use these images in this study. 

3.1.3. KP/512 



Observations by the KP/512 magnetograph ( |Livingston et al. 
1976) > were taken between 1974 and 1993. Full-disk images are 
2048 x 2048 pixels with a 1" pixel size, but with a generally sig- 
inificantly lower spatial resolution given by seeing. Continuum 
intensity images and magnetograms were built simultaneously 
from scans lasting approximately 40 minutes. Continuum inten- 
sity images are limited to a 4-bit intensity scale providing only 
15 discrete levels that make the separation of penumbra from 
umbra, or sunspots from limb-darkening, difficult. A brief de- 
scription of the procedure employed b y |Wenzler et al. ([2006 ) 



4.1 



Jones et al.| ( |1992| i and 



le magnetograms had a 



to resolve this is given in section 

enzler et al. (20061) both found that th 
background noise level of ~7-9 G; We nzler et al.| ( |2006| l found a 
mean of 8.1 G = cr noise 512 , which we adopt here. 

As the observatory is ground-based, seeing and weather con- 
ditions in general have a significant effect on the quality of the 
images. Work done by Wenzler et al.| ( |2006| l identified 1734 
dates with usable data over the period 1974 February 1 to 1992 
April 4. In 1992 there was a change of instrument and KP/512 
observations are only available as an overlap with the superced- 
ing instrument, KP/SPM, on 45 days between 1992 November 
28 and 1993 April 10. 



All modifications to produce level 1.8.2 magnetograms are de- 



|Wenzler et al.| (|2006 ) found that a conversion factor, fsi2, was 
necessary in order to multiply magnetic flux in KP/512 magne- 
tograms to bring them into agreement with KP/SPM. The official 
value of fsn quoted in the data release was 1.46, but Wenzler| 
et al. ( 2006[ l found that to get best agreement with PMOD a value 



of 1.60 was needed. The conversion factor allows other parame- 
ters, such as the magnetic flux level B t h, and B sat , to be the same 
in both KP/SPM and KP/512 reconstructions. We show, in sec- 
i 4.2 that fsn is dependent on B sat and has no freedom to vary 



tion 



independently. 



3.2. Observational TSI composite 

TSI observations are needed to fix the free parameter, B sat , in 
SATIRE-S. They also provide a basis with which to test the 
model reconstruction on various timescales. 

TSI observations began in 1978, but no single satellite 
has observed the whole period. Therefore, datasets have been 
combined into composites. There are currently three pub- 
lished composites adopting different approaches to combine 
the radiometric data and account for degradation, switch-offs 
and sudden glitches. These composites are known as PMOD 
(Frohlich 200 3), ACRIM ( jWillson & Mordvinov| [2003] l and 
IRMB ( |Dewitte et al.|2004[ ). In section|6|we compare the three- 
cycle TSI reconstruction with these composites. 

The current version of the PMOD composite, d41_62_1003 
2000) 1, covers the period 1978 to present. It is 
of data from four radiometers: NIMBUS7/HF, 
UARS/ACRIM II and the PM06V radiome- 



( Frohlich 



composed 
S MM/ AC RIM I, 
ter on the SoHO/VIRGO instrument (see references within 
[T998) . 



IFrohlich & Lean 



The ACRIM composite also uses 
NIMBUS 7/HF and~SMM/ ACRIM I, but where PMOD uses 
SoHO/VIRGO from 1996, the ACRIM composite contin- 
ues to use UARS/ACRIM II until 2000 before switching to 
ACRIMSAT/ ACRIM III. The latest version of the IRMB com- 
posite begins in 1984, employing NIMBUS7/HF, SMM/ ACRIM 
I, UARS/ACRIM II and then the DIARAD radiometer, on the 
SoHO/VIRGO instrument, from 1996. 

The current absolute irradiance level of P MOD is set at 
~1365 Wm" 2 during the solar minimum of 1986 (Crommelynck 
|et~aT1[T995] ). The latest TSI instrument, SORCE/TIM, has ob- 
served irradiance at ~4.5 WrrT 2 lower than VIRGO; recent cal- 
ibration measurements have confirmed this value to be more ac- 
curate (Kopp & Lean 2011). However, the relative variation of 
the PMOD composite is in good agreement with SORCE/TIM 
and it is the relative variation that is our main concern in this 
study; SATIRE-S can be scaled to either level. For these rea- 
sons, we normalise all three composites and SATIRE-S to the 
absolute level of TIM at solar minimum in 2008. 



3.3. SORCE/TIM 

The Total Irradiance Monitor (TIM), a radiometer on the 
SORCE satellite, has been observing TSI since February 2003 
with an absolute accuracy of 350 ppm, an instrumental noise 
of ~2 ppm and a stability of better than 10 ppm/yr (Kopp & 



Lawrence|[2005] |Kopp et al.||2005] |Kopp & Lean]|201 1| >. Apart 
from the lower absolute value of TSI, TIM also displays a shal- 
lower gradient than PMOD during the declining phase of cycle 



23. In this study we use version 11, level 3 six-hourly data^J av- 
eraged onto a daily cadence. 



scribed in more detail at http://soi.stanford.edU/magnetic/Levl.8/ 



SORCE/TIM data downloaded 201 1 July 1 1 
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By the time solar minimum had been reached, in December 
2008, TIM had been operating for about 5.5 years. This equates 
to a long-term l<x uncertainty in the change in TSI of ~55 ppm 
or ~0.07 WirT 2 . This is about the same as the difference in TSI 
variations seen between VIRGO and TIM over the same pe- 
riod and means that the gradients in TSI agree with each other 
within this uncertainty, not even taking into account the error on 
VIRGO (private communication, Greg Kopp). 



4. Homogenising KP/512, KP/SPM and MDI data 

In this study, we com bine KP/5 12 and KP/SP M w ith SoHO/MDI 
data. The studies of |Krivova et al.| (|2003 ) and |Wenzler et al.| 
(2004 2005 2006 ) also used these data, but since these publica- 
tions, MDI data have been recalibrated (see section |3". O) . This 
means that settings in SATIRE-S used by these authors may no 
longer be appropriate and need to be reconsidered. 

In this section, we discuss how the parameters in SATIRE-S 
are determined so that consistency is achieved across all three 
datasets. 



4.1. Sunspot areas and temperatures 



The study of Wenzler et al. 
KP/512 data. 



combined KP/SPM and 



They set sunspot thresholds so that agreement 
was found, in both datasets, with an independent ground-based 
composite record of sunspot area (Bal maceda et al.| 2005). The 
KP/SPM penumbral threshold was set in | Wenzler et al.| ( |2004| > 
at 0.92 so as to obtain good agreement with total sunspot area. 
|Wenzler| ( |2005| l confirmed that the ratio of umbral area to total 
sunspot area is 0.2 and used this to determine a KP/SPM umbral 
threshold of 0.70. 

In this study, we use the KP/512 sunspot data computed by 



Wenzler et al. (2006i for the KP/512 part of the reconstruction. 
The way these were obtained by |Wenzler et al.| ([2006 ) is as 
follows. Due to the 4-bit digitization of the KP/512 continuum 
intensity images, only the combined total umbral and penum- 
bral sunspot areas could be identified. Therefore, | Wenzler et al.| 
(2006) used the established ratio of umbra to sunspot area of 
0.2 to estimate umbral areas. Pores could not be identified, so 
that a simple multiplication of the total sunspot area brought 
the KP/512 sunspot areas into agreement with an independent 
sunspot area composite. 

Following from this, we vary MDI contrast thresholds to 
bring them into agreement with the spot ratio of 0.2 and 
KP/SPM. Using 847 continuum images from KP/SPM and MDI 
between 1999 and 2003, we compare areas and find that the 
penumbral and umbral thresholds should be 0.89 and 0.64 re- 
spectively. The difference between KP and MDI thresholds is 
partly due to the different wavelengths at which the instruments 
operated, as well as differences in stray light. 

The total effect of sunspots on TSI is degenerate between the 
temperature, T 4 , of the model atmosphere and the area of the 
identified sunspot. For consistency we use temperatures of 5450 
K for penumbra and 4500 K for umbra throughout, as in Wenzler 
|etaL| ( |2006> . 

4.2. The KP/512 magnetogram conversion factor, f 5n 



The introducion of f 51 2 (see section [3.1.3 i by Wenzler et al. 
( |2006| l to easily combine KP/512 and KP/SPM created an ad- 
ditional free parameter in the reconstruction. We are able to re- 



move the freedom in this parameter by considering how fsi2 re- 
lates to B sat . This is done as follows. 

Magnetic flux is registered differently in KP/512 and 
KP/SPM magnetograms, so we use the method of histogram 
equating ( Jones & Ceja|200T Wenzler et al.|2004 i to compare 
the signal observed in the magnetograms. In Fig. [1£, we com- 
pare unsigned magnetic flux in 18 KP/SPM and KP/512 mag- 
netograms taken on the same days in 1992/3. Crosses are the 
comparative magnetic flux calculated by sorting pixel values into 
ascending order and taking the mean in each of 1000 equally 
divided bins. By interpolating the KP/512 values onto a linear 
KP/SPM scale and averaging, the solid grey line results. This 
gives the corresponding mean level of magnetic flux in KP/SPM 
as observed by KP/512; the grey dashed lines show the standard 
error on the mean. Then, the mean level of magnetic flux over 10 
G intervals centred on each level of magnetic flux along the grey 
line is calculated. The ratio of these mean magnetic flux values 
in KP/512 to those of KP/SPM result in a conversion factor at 
every level of magnetic flux, as shown in Fig. [2^. In this plot, 
the solid line is the factor required to convert magnetic flux in 
KP/512 to bring it in to agreement with KP/SPM, as a function 
of magnetic flux. The dashed lines show the standard errors, also 
from Fig.[lJ. 

The official conversion factor value of f$\2 = 1.46 and the 
factor value of fsi2 = 1.60 found by | Wenzler et al.| ( |2006[ l are 
indicated by horizontal dashed lines in Fig. |2£. It now becomes 
clear where the official conversion factor and that of I Wenzlerl 



et al. ( 2006 1 come from: at high flux levels above 350 G the con- 
version factor is close to the official factor; below this it varies 
between 1.25 and 1.90, a range that very nearly averages to the 
value of |Wenzler et al.| ( j2006] l. 

In SATIRE-S, magnetic flux below B t h r is ignored and all 
magnetic flux above B sat is considered to have the same effect 
on irradiance, so only fluxes between these two values are im- 
portant in determining f^n. An approximate value for fsra can 
now be calculated by simply taking the mean of the conversion 
factor between B t hi-.5i2 and B sat 51 2. B t h ri 5i2 is fixed, therefore f^ 
depends only on the chosen value of B sati 5i2. The result is that 
fsi2 is completely determined by B satj 5i2, has ceased to be an in- 
dependent free parameter, and allows us to combine KP/SPM 
and KP/512 images while maintaining the same value of the sin- 
gle free parameter, B sat . 

Although B t hi-,5i2 = 21 G was determined using f^2 - 1.46, 
it changes very little with f 51 2, so we leave it at this value. In 
Fig. [2J5, the resulting conversion factor as a function of B sat> 5i2 
can be seen to lie in the range of B sat ,5i2 of 250 to 500 G (solid 
line) with the propagated standard errors (dashed lines) from the 
left plot at approximately + 0.04. In section 5.1 we find that the 
required value of B sat gi2 is 330 G. The resulting value of fsn, 



1.580, is close to the best value determined by Wenzler et al 



(2006) when the model was fitted to the PMOD composite. Now, 
however, fsi2 is determined independently of the PMOD com- 
posite. 

4.3. Removing dark high-flux magnetogram pixels 



Ball et al. ( |201 1\ found that high magnetic flux facular pixels 
identified in MDI magnetograms were incorrectly enhancing re- 
constructed irradiance relative to TSI observations. By ignoring 
pixels at and above this level, B cut MDi = 800 G, the agreement 
between modelled irradiance and SORCE/TIM improves, espe- 
cially on rotational timescales. Ball et al. ( |201 1\ incorrectly as- 
sumed that the physical reason was due to a decrease in con- 
tinuum contrast with increasing magnetic flux as in |Ortiz et aT[ 
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200 400 600 800 200 400 600 800 1000 

KP/SPM; Magnetic flux, G SoHO/MDI; Magnetic flux, G 



Fig. 1. Histogram equating plots, (a) 18 contemporaneous KP/SPM and KP/512 magnetograms and (b) 84 contemporaneous 
SoHO/MDI and KP/SPM magnetograms. The grey line is the interpolated mean of the black crosses with the dashed grey lines 
the standard error. Dotted lines highlight the relationship of B cut between the datasets. 
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Fig. 2. (a) The factor by which magnetic flux in KP/512 magnetograms should be multiplied in order to obtain the same magnetic 
flux as in KP/SPM magnetograms, as a function of magnetic flux (see text), (b) The best global conversion factor, fgi2, by which 
complete KP/512 magnetograms should be multiplied as a function of B sat 5i2, the SATIRE-S free paramter. Dashed lines in both 
panels represent the standard error as shown in Fig.[lJ. 



(2002). According to Kok Leng Yeo (private communication), 
the magnetic flux of pore and sunspot canopies are being de- 
tected near the limb and falsely attributed as faculae. The im- 
proving effect of removing these pixels appears to be most sig- 
nificant at ~550 G, but down to this level real magnetic features 
away from the limb are also removed. Therefore, the more con- 
servative cutoff level of 800 G is maintained in this paper. 

For consistency, this result needs to be propagated through to 
KP/SPM and KP/512 magnetograms. Figure [TJi is constructed 
in the same way as Fig. [T^, as described in the previous sub- 
section. To calculate the shift, B cuti spM, we compare 84 magne- 
tograms observed within two hours of corresponding MDI mag- 
netograms between 1999 and 2003. Using such close temporal 
proximity reduces the effect of evolving active regions. This in- 



dicates that when B cut mdi = 800 G then B cut; spM = 395 G. We 
repeat the same process for KP/512 using Fig.[T^. It follows that 
B cu t,5i2 = 616/f5i2 G, where fsn is the conversion factor (see sec- 
tion [3X3). 

We have established that B cllt is dependent on the chosen 
value in MDI. While this was a somewhat arbitrary choice, it 
is not a freely varying parameter. Since the freedom of f^n can 
be removed and made dependent on B sat (see section 4.2 1, no 
additional free parameters have been introduced here. 



4.4. Correcting pre-1990 KP/512 magnetograms 



Arge et al.| ( |2002| l reported that, prior to 1990, the total photo- 
spheric magnetic flux from KP/512 Carrington maps (NSO) was 
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Fig. 3. (a) Total photospheric flux, obtained from KP/512 Carrington maps, before applying a correction (black dots), after either 
a scaling (green) or a shift is applied (blue), (b) SATIRE-S reconstruction of TSI using KP/512 continuum intensity images and 
magnetograms only, with either a scaling or a shift compared to the ACRIM and PMOD composites. All datasets are smoothed and 
normalised to the minimum of 1986. Gaps in the curves are present when data gaps are larger than 27 days. See main text for scaling 
and shifts applied to KP/512 Carrington maps and magnetograms. 



lower compared to that from Wilcox (WSO) and Mount Wilson 
(MWO) Solar Observatories Carrington maps. They used their 
analysis to develop a correction procedure, but this was only ap- 
plied to a limited number of magnetograms, most of which were 



recorded during inactive times of the solar cycle. Wenzler et al 



(2006) needed to use all available KP/512 magnetograms. They 
found a straight-line fit between the corrected and uncorrected 
total Carrington-rotation integrated photospheric flux, (see Fig. 
3 of [Wenzler et al.[|2006| ). The linear fit was excellent (corre- 
lation coefficient, r c =0.98), so all pre- 1990 magnetograms were 
scaled by the regression slope value of 1 .242 (the intercept of the 
fit was 0.306). 

In Fig. [3^ we show the effect of scaling by a factor of 1 .242 
(green dots) compared to uncorrected Carrington maps (black 
dots), shown as the total photospheric flux of each Carrington 
roration between 1974 and 1992. The uncorrected Carrington 
maps can be seen compared to the WSO and MWO in Fig. 1 
of |Arge et aLlp002] l. Figure 5 of | Wenzler eTaT] ( |2006> shows 
the effect after scaling KP magnetograms by a factor of 1.242 
and while the scaling seems to have improved the NSO results 
during the 1974-1976 minimum relative to WSO and MWO, it 
appears that the effect is too large between 1987 and 1989. The 
effect on the maximum of cycle 21 seems to be reasonable but in 
general the scaled level sits above MWO and WSO level of total 
flux at all other times. This may be due to the fact that the scaling 
factor is mainly derived from Carrington maps recorded during 
solar activity minima. Extrapolating the scaling factor to more 
active periods may cause magnetic flux at cycle maximum to be 
overestimated. Therefore, we also consider an alternate correc- 
tion to apply. 

While a general correction to pre- 1990 data would likely take 
the form of a time-dependent shift and multiplication, this would 
require a detailed analysis of all magnetograms. In this paper, 



we consider the effect of applying both a constant shift or a scal- 
ing. To calculate, we take the difference between uncorrected 
and corrected total NSO flux below 3 xlO 15 Weber during low 
activity times. The value of 3 xlO 15 Weber is arbitrary but en- 
compasses the majority of low activity Carrington maps. The 
mean difference is 4.7 G xlO 14 Weber. The effect of applying 
this shift to the Carrington maps is also shown in Fig. [3^ as blue 
dots. At solar minima, the shift results in similar values to the 
scaling, but the cycle amplitude is lower and agrees better with 
MWO and WSO during the early part of cycle 22 maximum. 
Unfortunately, the shift derived for the synoptic charts cannot 
be converted to a shift in the KP/512 magnetograms: total mag- 
netic flux reported in KP/512 magnetograms can exceed NSO 
Carrington maps by a factor of well over 100 due to massive 
flux cancellation in the lower-resolution Carrington maps. 

Instead, we take a different approach to calculate a con- 
stant shift for the pre- 1990 KP/512 magnetograms. First, we as- 
sume that the scaling factor of 1.242 found by Wenzler et al. 
( 2006 ) is appropriate for low activity magnetograms. We can use 
the difference between low-level magnetic flux in uncorrected 
and scaled magnetograms as a constant shift. For this we use 
the threshold level determined by Wenzler et al. (2006) for the 
KP/512 magnetograms, i.e. B t hr,5i2 = 3cr ma ^ 51 2 = 3 X 8.1 G = 
24.3 G. This value was found by using only magnetograms dur- 
ing the solar minimum of 1986 and, therefore, is appropriate for 
use in calculating a constant shift. The difference between scaled 
and unsealed Bthr,5i2 is 5.9 G. In Fig. [3J?, reconstructions using 
KP/512 magnetograms only are shown, with scaled in green and 
shifted in blue. All datasets in this plot are shown as smoothed 
time series and normalised to the solar minimum of 1986. All 
smoothings of plotted time series throughout the whole paper, 
and for detrending data, are performed with a Gaussian win- 
dow with an equivalent boxcar width of 135 days. The result 
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Fig. 4. Mean background magnetic flux distribution of magnetograms from 1996-7 (left) and 2008 (right). Both plots have the same 
colour scale and contours are in 1 G levels. 



of adding a constant instead of scaling is to reduce the cycle am- 
plitude prior to 1990. 



considering the range of uncertainty on the reconstruction, we 



include both the shift and scaling as discussed in section 5.2 



It is not immediately clear which approach, scaling or shift- 
ing, should be adopted. Therefore, we use the TSI composites 
as a guide, assuming that when they agree they are a reasonable 
estimate of the true value of TSI at that time. In Fig. [3j), the 
ACRIM (yellow) and PMOD (red) composites are also shown. 
The IRMB composite is no longer extended back to 1978 and 
so is not considered here. For the period of time covered in this 
plot, ACRIM and PMOD use the same TSI data sources post- 
1983, though the way they have been combined is different and 
after 1989, the trends clearly differ. Therefore, the agreement 
over 1984 - 1989 is expected. Even though ACRIM and PMOD 
use different datasets prior to 1984, the two composites agree 
well between 1981 and 1983. Before 1981, instrument degra- 
dation corrections are applied to PMOD, but not to ACRIM. 
Neither scaling nor shifting fully corrects the KP/512 magne- 
tograms (private communication, Jack Harvey), though together 
they should give a good indication of the expected range of plau- 
sible reconstructions, with shifting minimising the cycle ampli- 
tude and scaling maximising the cycle amplitude. As the recon- 
struction with a shift applied clearly agrees better with ACRIM 
and PMOD over the 1981 to 1989 period than using the scal- 
ing, we choose to adopt the shift for the final reconstruction (see 
section [§}. This, then, assumes that shifting is a more appropri- 
ate approach to reconstructing TSI than using a scaling. When 



4.5. SoHO/MDI magnetograms 

The two KP datasets, KP/512 and KP/SPM, have been set up 
to use the same parameters. But to combine MDI with these 
datasets, an overlap period is needed to find the appropriate B sat 
through linear regression. It is, therefore, important that the im- 
ages used in SATIRE-S are stable over time. For MDI this does 
not appear to be the case. 

The lowest level of magnetic flux in the MDI magnetograms, 
Bthr.MDi = 3cr noise MDI is dependent on disk position. The changes 
to the MDI magnetograms from level 1.5 to 1.8.2 required 
a re-examination of the mean background flux to recalculate 
o"noise,MDi as a function of disk-position. 

d2002) and ap- 



We used a similar approach tolOrtiz et al. 



plied this to two sets of 29 magnetograms from both 1996/7 and 
2008, chosen for the minimal activity present on the disk. The 
procedure is as follows. A 100 x 100 pixel box centred on every 
pixel in the magnetogram was used to determine if the value of 
the central pixel exceeded three standard deviations of the mean. 
If the pixel did exceed this threshold, it was set to the mean. 
This process was repeated two further times and then the me- 
dian of the 29 distributions was taken. The contour and surface 
plots of the level 1.8.2 magnetograms in Fig.|4]show the mean 
background flux distribution in 1996-7 (left) and 2008 (right). 
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Two things are clear from Fig.|4j the distribution is not uni- 
form across the disk and it changed between 1996 and 2008. 
The non-uniformity across the disk, seen as an increase from the 
upper-left to lower-right in both plots, can possibly be explained 
by the inhomogenous nature of the Michelson interferometers 
in the field of view: slightly different wavelengths are observed 
across the disk. However, the increase over time of the lower- 
right quadrant is of most concern here. Let us now consider pos- 
sible causes for such a change. 

Aside from several shorter periods without observation, there 
were two extended periods, between June and October 1998 and 
from December 1998 to February 1999, when SoHO lost con- 
tact and was in an uncontrolled state. During the early years of 
operation, prior to this loss of contact, the focus setting of MDI 
was deliberately set to ensure the quality of the 'high resolu- 
tion field' mode observations and, as a result, 'full-disk' mode 
was slightly defocused. When contact was re-established, fol- 
lowing an analysis that showed this defocusing was not required, 
the instrument was fully refocused, which significantly changes 
the point spread function (PSF). Subsequent periodic refocusing 
was performed to correct a drift and maintain the focus, which 
should basically maintain the PSF at roughly the level achieved 
after the loss of contact. This, among other effects such as e.g. 
the front window darkening, thermal effects, a small change in 
the observed wavelength and any unknown effects during the 
period without contact, may have contributed to a change in the 
registered magnetic flux as a function of disk-position (private 
communication, Rock Bush). We found that dates when focus 
changes were performed correlate with a change in distribution 
of flux in the lower right quadrant. Indeed, by bringing the in- 
strument into focus and narrowing the PSF, the expected out- 
come would be an increase in the observed magnetic flux, and 
this is exactly what can be seen in Fig. [4] The changes warrant a 
deeper investigation of the magnetogram changes (see forthcom- 
ing thesis by Will Ball). Here, we are mainly interested in global 
sensitivity to the magnetic flux through the effect they have on 
a reconstructed irradiance, which can be judged by comparing 
reconstructions using MDI and KP/SPM data prior to and after 
the loss of telemetry. 

In Fig. [5] daily and smoothed reconstructions for the overlap 
period of MDI (blue) and KP/SPM (orange) are shown. The two 
reconstructions have a unity regression (see section |5.1| l post- 
SoHO gaps (grey shading), but show a shift prior to the gaps. It 
can be rectified by a simple shift in the MDI-based recnstruction 
of 0. 14 Wm 4 . The shifted reconstruction is shown as daily (thin, 
black) and smoothed (dashed) lines. The agreement of the trends 
in the shifted MDI and KP/SPM reconstructions is excellent. 

Note that a defocus also reduces contrast in the continuum 
intensity images (Danil ovic et al.||2008) , which would also in- 
troduce a bias between sunspot and umbral areas if determined 
from the period prior to loss of contact relative to determining 
it afterward. We chose thresholds based on continuum intensity 
images from after the loss of telemetry (see section 4. 1 1, so this 
is not a concern here. 

With both of these results, along with the knowledge that the 
source of the difference almost certainly lies with MDI given the 
known changes over the SoHO data gaps, we decide to only use 
MDI data from after the data gaps and use KP/SPM to recon- 
struct prior to this. 

5. Reconstruction of TSI for cycles 21-23 

In this section, we explain how we produce the three-cycle re- 
construction, independently of any TSI composite after 1989, 



and how the errors have been formulated. The ultimate aim here 
is to produce a composite-independent reconstruction of TSI. 
The SORCE/TIM data, unused in the composites at the time of 
writing, is therefore an ideal candidate for use as an independent 
basis to fix the model's free parameter. 

5.1. Combining MDI and KP/SPM data 

SATIRE-Smdi is normalised to TIM for 2105 days between 2003 
February 25 and 2009 October 31. B sat MDi is varied until the 
best fit regression with a slope of 1.00 is found at 443 G, using 
the FITEXY linear regression routine (Pre ss et al.| [T992). This 
is shown as the solid line in the left regression plot of Fig. [6] 
The correlation coefficient at this value of B sat mdi is r c tim mdi = 
0.979. 

Note that a 30-day period centred on the huge sunspot group 
passage that reached disk centre around 2003 October 29 has 
been omitted from this fit. We find that the fit is relatively sta- 
ble with B sat MDi = 443 G for the entire overlap period with TIM 
except during this passage. SATIRE-S generally reproduces the 
influence of sunspots well, but this sunspot group seems to have 
been so exceptionally large that the assumed sunspot properties 
appear to have changed. SATIRE-S significantly overestimates 
the darkening produced by this sunspot group and this leads to a 
change in B sat mdi that causes the rest of the time series to be less 
well reproduced. The spot is only removed when determining 
B sat ,MDi, but considered in the comparisons with the TSI com- 
posites (section[6]). 

Next, we combine SATIRE-Sspm with SATIRE-Smdi using 
the overlap period of 895 days between 1999 February 18 and 
2003 September 21. SATIRE-Sspm is normalised to SATIRE- 
Smdi- Using SATIRE-Smdi with B sat ,MDi = 443 G, a unity re- 
gression fit for SATIRE-Sspm is achieved when B sat s PM = 330 
G. The regression plot is shown on the right-hand side of Fig. [6] 
The correlation coefficient of this fit, r c mdlspm, is 0.943. 

With SATIRE-Sspm established, the link to SATIRE-S512 is 
automatically prescribed. The same normalisation factor found 
to bring SATIRE-Sspm in to line with SATIRE-Smdi is used for 
SATIRE-S512- The conversion factor, fsn, for B satj 5i2 = 330 G 
(see Fig.[2]3) is 1.580, to the nearest 0.005. The entire time se- 
ries is normalised to TIM at the December 2008 minimum as a 
reference period. 

Note that the method to combine MDI with KP/SPM differs 
to that of combining KP/SPM with KP/512. In the latter case, 
magnetograms from both datasets are compared and adjusted 
so that the same value of the free parameter can be used. This 
method is not possible with MDI because the spatial response is 
not uniform, as discussed in section [43] 

We use the TSI composites to justify the addition of a con- 
stant to pre-1990 KP/512 magnetograms (see section 4.4 1 and, 
therefore, the reconstruction is independent of the composites 
after 1989. We denote this full reconstruction as SATIRE-Si n d. 
In Fig. [7] the complete SATIRE-Si n d reconstruction is shown as 
daily (thin lines) and smoothed (grey lines) time series. All that 
remains is to estimate an uncertainty on this result. 



5.2. Uncertainty estimation on the reconstruction 

While errors associated with the continuum intensity images 
and magnetograms enter into the SATIRE-S reconstruction, the 
largest uncertainties arise in the calculation of the emergent in- 
tensities for the different components, in particular for the facu- 
lae. These uncertainties are difficult to quantify, though it turns 
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Fig. 5. Reconstructions of KP/SPM (orange) and SoHO/MDI (blue) for the overlap period May 1996 to September 2003; KP/SPM 
is normalised to MDI after February 1999. Daily data are shown as thin lines, smoothed as thick, with the grey boxes covering 
periods when SoHO was out of contact. Prior to the loss of contact, MDI is also plotted shifted by a constant value (daily data, black 
thin solid line; smoothed, dashed thick line). 



1363 
1362 

1361 



CO 

I 

Ld 
0£ 

I— 
< 

in 



1358 



_ I 

: r c,TIM:MDI — 


11 11 I 

0.979 


+/ 


L $,TIM:MDI = 


2105 








: + 






: +, 


+ 




: A- 






/ 







SORCE/TIM, Wm" 



1365 r 



1358 1359 1360 1361 1363 

-2 



tn 
CO 
I 



< 
CO 



r c,MDI:SPM — 0.943 



1 364 - ff,MDI:SPM 

1363 ~r 



895 



1362 - 
1361 - 
1360 ~r 



1359 V 

1359 1360 1362 1363 1365 
SATIRE — S MDh Wm" 2 



Fig. 6. Regression of (left) TIM and MDI and (right) MDI and KP/SPM. The regression line is shown by the solid grey line and lcr 
uncertainty is represented by the dashed grey lines. 



out that changes in the facular contrast are largely absorbed by 
changes in B sat as far as TSI reconstructions are concerned. 

We therefore try and estimate the uncertainty in the SATIRE- 
S reconstructions by exploring the change in the modelled TSI 
that results from variations in the free parameter, B sat , that fall 
within the allowed fit parameters, such as e.g., quantified by a 
lcr change in the regressions. We use the difference between the 
datasets being regressed as an estimate of the standard deviation 
for both datasets in the FITEXY routine. 



Starting with SATIRE-Smdi we found that, the formal lcr 
error on the unity fit gradient with SORCE/TIM, ±0.006, corre- 
sponds to ±5 G in B satMDI . This uncertainty also encompases the 
best value for B sat mdi when a least-squares fit, rather than a re- 
gression, is used to find agreement between TIM and SATIRE-S. 

The uncertainty over the KP/SPM part of SATIRE-S Ind is 
determined from the uncertainty range on the unity regression fit 
between SATIRE-Smdi and TIM. In other words, SATIRE-S SPM 
is fitted to SATIRE-Smdi using B sat , M Di = 438 and 448 G for the 
upper and lower boundaries, respectively. Then, the lcr errors 
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on these unity regression fits are used to find an expanded uncer- 
tainty range prior to 1999. The lcr error on the unity regression 
gradient is 0.017 and the resulting B sat ,spM is 307 G and 365 G 
for the upper and lower range, respectively. 

The uncertainty range of B sat> 5i2 in SATIRE-S512 includes the 
standard error of approximately +0.04 on the associated conver- 
sion factor, fsi2 (see section |4~2] >. We take values of fsn that ex- 
acerbate the uncertainty range. For B sat ,5i2 = 365 G, this equates 
to f 5 i2 = 1.600 and for B satj5 i2 = 307 G^ f 5 i 2 = 1.550. 

Finally, we must consider the effect of using a scaling or con- 
stant shift correction to KP/512 magnetograms prior to 1990. We 
consider both approaches so that on each date the pre- 1990 cor- 
rection that produces the largest deviation from SATIRE-Si n d, 
either highest or lowest in value, is adopted as the upper and 
lower uncertainty, respectively, to maximise the range. 

The solar minimum of December 2008 is used as the ref- 
erence date, assuming no long-term uncertainty at that time. 
Therefore, the two extreme reconstructions that constitute the 
uncertainty range are normalised to that time, though this does 
not mean that there is no uncertainty on the absolute value of TSI 
then. The resulting reconstruction is shown in Figj7] along with 
the smoothed uncertainty range shown as a contour around the 
smoothed dataset. While SATIRE-Si„d pre- 1990 is not indepen- 
dent of the TSI composites, the uncertainty range is as the uncer- 
tainty includes both the maximum range produced by shifting or 
scaling KP/512 magnetograms prior to 1990. 

6. Comparison of SATIRE-Si n d with TSI Composites 

In this section we compare SATIRE-Si n d to the three TSI com- 
posites in detail. This returns us to a debate over which compos- 
ite (ACRIM, IRMB o r PMOD) provides the best record of TSI. 
Wenzler et"aL] ( |2009| l used S ATIRE-S to do this by optimising 
the model for each composite and found that SATIRE-S fitted 
to PMOD gave the best correlation coefficients, while ACRIM 
and IRMB displayed trends that PMOD and SATIRE-S did not 
show. 

SATIRE-Si n d provides a much stronger test as it covers 
a longer period and is independently optimised. In Fig. [8] 
all data are shown together, normalised to the mean value of 
SORCE/TIM (not shown) over the six month period centred 
on the December 2008 solar minimum. The data have been 
smoothed to aid clarity, with PMOD in red, IRMB in green, 
ACRIM in yellow and SATIRE-S Ind in blue. In Fig. [9) the dif- 
ference between the composites and SATIRE-Si n d are shown 
for daily data for the period 1978-2009. Before 1984, data are 
shown on each date that data are available in PMOD, ACRIM 
and SATIRE-Si n d; after 1984, data are shown when available in 
all four datasets. Table [T] summarises comparisons between the 
TSI composites and SATIRE-Sind- It is divided in two halves 
deliniated by the double line: the top half shows the results of 
using the daily dataset while the bottom half is for detrended 
data only so that the rotational variability alone can be consid- 
ered. Detrending in all cases has been done by subtracting the 
smoothed dataset from the original. These two halves are further 
divided into comparisons spanning 1978-2009 and 1984-2009, 
the latter period being the one to which the IRMB composite is 
restricted. All comparisons use only dates when data are avail- 
able in all datasets. This means that the correlation coefficients, 
r c , and slopes listed in Table [T] are strongly weighted by cycle 
23, as most data points in SATIRE-Si„d are during that period. 

In Fig. [8} PMOD displays a smaller trend than the other two 
composites at nearly all times. Also, SATIRE-Si n d is in good 
agreement with PMOD at all times, except around the peak of 



cycle 22. SATIRE-Si„ d agrees better with PMOD than with other 
TSI composites and better than the composites agree with each 
other, reflected in the 1984-2009 daily non-detrended results. 
Since the end of cycle 21, SATIRE-Sind shows a decrease in 
minima levels within the uncertainty range. This secular change 
appears to be most significant between the minima of cycle 23, 
in agreement with PMOD, and does not show the more marked 
secular variations as in the IRMB and ACRIM composites. 
Interestingly, IRMB and SATIRE-Sind appear to have the best 
agreement during the declining phase of cycle 23, with a gra- 
dient in better agreement with SORCE/TIM than either PMOD 
or ACRIM (for SATIRE-S this is by construction). Note that the 
regressions and correlation coefficients improve in the shorter, 
later period from 1984, consistent with the fact that the older 
data prior to 1984 are less certain. 

Considering the detrended TSI comparisons between 1984 
and 2009, which are measures of rotational modulation of ir- 
radiance and the effects of active region evolution, agreement 
between all datasets is generally good. IRMB has the best 
agreement with SATIRE-Sind, followed by ACRIM and finally 
PMOD. The agreement between SATIRE-Si nd and IRMB is al- 
most as good as the best agreement between TSI composites, 
which is found between IRMB and ACRIM. It is also clear that 
the rotational scatter in Fig. [9] between IRMB and SATIRE-S^d 
during cycle 23 is much lower than with the other TSI compos- 
ites though the longer-term trend during the rising phase is not 
the same, especially before 1997; this entire period utilises the 
DIARAD instrument on SoHO. 

The high level of agreement in the detrended time series 
between SATIRE-Si„ d and IRMB means that SATIRE-Sind is 
able to recreate rotational variability to the level of uncertainty 
between the TSI composites. This goes some way to explain- 
ing why IRMB has a relatively good regression gradient with 
SATIRE-Sind when considering the non-detrended data, despite 
the large disagreement in the smoothed long-term trend shown 
in Fig. |8] 

The regression slopes also indicate that the rotational varia- 
tion in SATIRE-Sind falls in between the TSI composites, with 
ACRIM having the largest amplitude, and PMOD and IRMB 
the lowest. It should be borne in mind that these results may 
be affected strongly by active regions. It is well documented 
flUnruh et alT| |2008 ; Bal l et al.pOTT) that TSI can be overesti- 
mated in SATIRE-S when large facular regions transit the edge 
of the solar disk. Some of this overestimation has been im- 
proved upon through the modification described in section |4.3| 
It may be the case, though, that SATIRE-S^d still overestimates 
rotational varability and that is why the regression slopes sug- 
gest that SATIRE-Si n d has a higher amplitude than PMOD and 
IRMB. 

In general, the composites and reconstruction agree well 
with each other on rotational time scales. The major cause of 
differences between composite datasets are likely to be produced 
by the use of different data sources in cycle 23. Assuming, for 
the sake of argument, that SATIRE-S^d is a good representation 
of TSI behaviour we find that PMOD is the most reliable com- 
posite of TSI, taking into account rotational, cyclical and secular 
variations, although on rotational time scales IRMB is superior. 

Given that PMOD is the composite by far the closest 
to SATIRE-Sind overall, we concentrate in the following on 
an in-depth comparison of the twoS.S.Dhomse@leeds.ac.uk. 
Figure [10] shows in the top panel PMOD (light red, daily; red, 
smoothed) and SATIRE-Si„d (light blue, daily; blue, smoothed) 
irradiance. The bottom plot shows the difference between the 
two, PMOD minus SATIRE-Sind (black, smoothed). Also in- 
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Fig. 7. The entire SATIRE-Si n[ j dataset. Daily data are shown as solid grey lines when dates are contiguous. The thick black line is 
smoothed data and the uncertainty range is shown, smoothed only, as thin black lines. Gaps in the curves are present when data gaps 
are larger than 27 days. 

Table 1. Comparison of SATIRE-Si n d with TSI composites using the original time series (upper half of the table) and detrended 
(lower half). The table is further split into two to allow a fair comparison with the shorter IRMB composite. Comparisons are made 
only for dates when data in all datasets are available. 





Daily dataset comparison (1978 


- 2009) 




Series 1 


Series 2 




Slope 


No. points 


ACRIM 


SATIRE-Si n d 


0.879 [0.772] 


0.832 


5663 


PMOD 


SATIRE-Si n d 


0.936 [0.875] 


1.000 


5663 


ACRIM 


PMOD 


0.866 [0.750] 


0.830 


5663 




Daily dataset comparison (1984 


- 2009) 




ACRIM 


SATIRE-Si n d 


0.903 [0.815] 


0.835 


5150 


IRMB 


SATIRE-Si n d 


0.924 [0.854] 


0.978 


5150 


PMOD 


SATIRE-Si„d 


0.942 [0.888] 


0.992 


5150 


ACRIM 


IRMB 


0.924 [0.853] 


0.856 


5150 


PMOD 


IRMB 


0.922 [0.850] 


1.013 


5150 


ACRIM 


PMOD 


0.894 [0.799] 


0.841 


5150 


Daily dataset comparison, detrended (1978 - 2009) 


ACRIM 


SATIRE-Si n d 


0.866 [0.750] 


0.986 


5663 


PMOD 


SATIRE-Si„d 


0.873 [0.762] 


1.039 


5663 


ACRIM 


PMOD 


0.830 [0.689] 


0.947 


5663 


Daily dataset comparison, detrended (1984 - 2009) 


ACRIM 


SATIRE-Smd 


0.889 [0.791] 


0.992 


5150 


IRMB 


SATIRE-Si„d 


0.940 [0.883] 


1.036 


5150 


PMOD 


SATIRE-Si„d 


0.877 [0.769] 


1.032 


5150 


ACRIM 


IRMB 


0.945 [0.892] 


0.957 


5150 


PMOD 


IRMB 


0.911 [0.831] 


0.994 


5150 


ACRIM 


PMOD 


0.847 [0.718] 


0.960 


5150 



dicated are the smoothed uncertainty range for SATIRE-Si n d, 
as thin, light blue and black lines in top and bottom plots re- 
spectively, and PMOD black and red error bars from Frohlich 
(2009). The values of these errors are given in the middle table 
of Table [2] These error bars were derived by Frohlich (2009} as 
an estimate of the long-term stability in the instruments and the 
calculation is described within that publication. The error bar at 
the 1986 minimum is very small because this minimum was used 
as the reference date for comparison with other minima. In the 
plot, PMOD is normalised to the mean of SATiRE-Si n( j at the 
2008 minimum, 2008 December 15 +90 days, to help discuss 



the change between cycle minima and bring irradiance levels in 
line with SORCE/TIM. Dashed lines in blue and red guide the 
eye to minima levels in 1986, 1996 and 2008. 



The immediate impression is that across the 5899 days, 
where data are available in both time series, there is very good 
agreement on all timescales for the majority of the time period. 
The long-term trend is very similar over the 3 1 year period of 
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1975 1980 1985 1990 1995 2000 2005 

Dote (yeor) 



Fig. 8. Comparison of smoothed TSI composites and SATiRE-Si n d, all normalised to SORCE/TIM at the minimum of December 
2008. ACRIM is in yellow, IRMB in green, PMOD in red and SATIRE-Si n[ j in blue. Gaps in the curves are present when data gaps 
are larger than 27 days. 




1980 



1985 



1990 



1995 
Dote (year) 



2000 



2005 



Fig. 9. Difference plots of (top) ACRIM - SATIRE-S Ind , (middle) IRMB - SATIRE-S Ind , (bottom) PMOD - SATIRE-S Ind . The 
plots cover the period 1984-2009 for dates when data are available in all four datasets and prior to this when data are available in all 
datasets excluding IRMB. 



overlap and the comparison shows that SATIRE-Si n d can explain 
almost 92% of variability in the TSI composite]^] 

In the lower part of Table[2] comparisons have been made for 
the full cycle and each individual cycle. The agreement is high- 
est in cycle 23 during which S ATIRE-Si n d is able to explain more 
than 96% of the variation. We note that the difference plot shows 
that the gradient in TSI during the rising and falling phases of cy- 
cle 23 is slightly lower in SATiRE-Sr n d, a consequence of using 
TIM to fit the reconstruction, which shows a gentler decline of 
TSI with time than PMOD. 

A significant result is that SATIRE-Si n d is able to reproduce 
the secular, inter-cycle, change in PMOD TSI, mimicking the 
-0.2 WtrT 2 decline between 1996 and 2008. There is a slight dif- 
ference between minima values in 1996, relative to 2008: within 
the uncertainty range, SATiRE-Si n d suggests a decline between 

cycle 23 minima of between 0. 1 1 and 0.32 WtrT 2 , encompass- lower panel of Fig. [TO] This is likely to be the result of the poorer 
ing the result from PMOD. The best fit from SATIRE-Sind also data quality of KP/512 and the larger uncertainty in PMOD dur- 
ing the first two cycles. The increased difference in 1990/1 oc- 
curs in KP/512 data in the period where no additive correction is 
applied to the magnetograms (see section 4.4 1. There is an impli- 
cation that a scaling correction per se should not only be applied 



falls within the error bars of PMOD in 1996. SATIRE-Si nd agrees 
with PMOD on the secular change over cycle 22 and possibly 
shows a slight secular increase over cycle 21 as can be seen in 
Table [2] though data are sparse prior to 1977 and there is no 
change within the uncertainty range. These results, summarised 
in Table [2] calculate the 1976 minimum using the mean of data 
between June 1975 and July 1976, while all other maxima and 
minima are for six month periods centred on the maximum and 
minimum months of the sunspot cycle from the Solar Influences 
Data Centre (SIDC) smoothed monthly sunspot dataset ( |SIDC-| 
|team|2011| >. 

Although the long-term agreement between PMOD and 
SATIRE-Si n d is generally very good, there are periods of one to 
two years during the maxima of cycles 21 and 22 where the dif- 
ference between them reaches ~0.3 WtrT 2 , as can be seen in the 



3 Note that the correlation coefficients and regression slopes between 
PMOD and SATIRE-Si nd differ in Tables[T]and^]due to a change in the 
data points considered. 
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Table 2. Detailed comparison between PMOD and SATIRE-S considering only dates when both datasets overlap 1 . 



Cycle Maximum 2 and Minimum TSI, Wm 2 


Series 


MinC20/21 J 


MaxC21 Min C2 1/22 Max C22 Min C22/23 


Max C23 


Min C23/24 


SATIRE-Si„ d 


1361.06 


1362.01 1361.08 1361.96 1361.01 


1361.74 


1360.82 


PMOD 




1362.02 1361.14 1362.10 1361.07 


1361.82 


1360.82 


Cycle Maximum and Minimum TSI Relative to 2008 Minimum, Wm 1 


Series 


Min C20/21 


MaxC21 Min C2 1/22 Max C22 Min C22/23 


Max C23 


Min C23/24 


SATIRE-S^ 
PMOD 


24 +u 14 

"•^-0.06 


, J9+0.4U 26 +U14 1 14 +042 20 +(m 

L - Ly -Q.M "■■'"-0.06 1- -0.09 "■•'"-0.09 

1.21 0.32±0.01 1.28 0.25±0.10 


92 +0 02 

"• 7 -0.02 

1.00 


00 +o ou 
"■""-o.oo 

0.00±0.16 



Statistics of Full and Cycle Periods 



Period 


Start date 


End date 


r c [r;.] 


Slope 


No. points 


Full 


1978/12/10 


2009/10/31 


0.958 [0.919] 


0.995 


5899 


Cycle 23 


1996/05/15 


2008/12/15 


0.981 [0.963] 


0.978 


3760 


Cycle 22 


1986/09/15 


1996/05/14 


0.917 [0.840] 


1.040 


1148 


Cycle 21 (part) 


1978/12/09 


1986/09/14 


0.898 [0.806] 


1.050 


740 



Notes. 1 The uppermost table shows the absolute values at activity cycle maximum and minimum. The middle table gives irradiance values 
relative to the December 2008 minimum with uncertainty estimates from SATIRE-Si„d and errors from PMOD. The bottom set give the regression 
slopes and correlation coefficients from the different cycles of the reconstruction. 2 TSI values are calculated from ±90 days centred on the dates of 
maximum and minimum smoothed monthly sunspot number from SIDC data. 3 Only SATIRE-Si n d data exist around the minimum of cycle 20/21 
minimum. Minimum values are estimated between June 1975 and July 1976 due to the lack of data around the sunspot minimum of June 1976. 



PMOD doily 
PMOD smoothed 
SATIRE-S,„ a doily 
5ATIRE-S lna smoothed 
SATIRE-S,„ a uncert. 




SATIRE-S lnJ daily 
PMOD - SATIRE-S|„ d smoothed 
PMOD - SATIRE-S.„ a uncert. 



1980 



1985 



1990 



1995 
Date (year) 



2000 



2005 



Fig. 10. In the upper plot, PMOD (daily data, light red; smoothed, thick red line) and SATIRE-Si nf j (daily data, light blue; smoothed, 
thick blue) between 1978 and 2009 normalised to SORCE/TIM at December 2008. The thin blue lines mark the uncertainty range 
of SATIRE-Si n d (only smoothed values plotted). In the lower plot, the difference between PMOD and SATIRE-Si n( j is shown (daily, 
grey; smoothed, black) along with the difference of the uncertainty with respect to PMOD. The black and red error bars are the 
errors from Frohlich (2009) in the upper and lower plots, respectively. Dotted vertical lines indicate cycle maxima and minima. 
Dashed horizontal lines signifying cycle minima are plotted to aid the reader. Gaps in the curves are present when data gaps are 
larger than 27 days. 



to pre-1990 magnetograms. Only applying a constant shift as a 
correction is probably also not correct and a combination of the 
two is probably more appropriate. A fuller investigation should 
be performed to understand the time-dependent problems in the 
KP/512 magnetograms, expanding on the work done by Arge 
|etaL| ( |2002> . 

There is also clearly a different six-monthly trend appar- 
ent in the period from late-1996 to 1998 where the gradient 



in SATIRE- Si n d is lower than in PMOD. Although MDI data 
over this period have not been employed for the final reconstruc- 
tion, for this small stretch both SATIRE- S M di and SATIRE-S S pm 
show a very similar trend (see Fig. BJ. |Krivova et al.| ( |2011a| l 
discuss this difference in some detail and also compare other 
radiometric observations that show significant deviations from 
PMOD. It is therefore very likely that the degradation correc- 
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tions applied during this period to VIRGO data may profit from 
being revisited. 

7. Discussion and conclusions 

In this paper we present the first reconstruction using the 
SATIRE-S model to fully cover solar cycles 21, 22 and 23. 
The reconstruction is compared with three TSI composites con- 
struc ted by|Willson & Mordvinov| ( |2003| l, |Dewitte et aT[ ([2004 ) 
and |Frohlich & Lean| ( |1998| l and referred to as ACRIM, IRMB 



and PMOD respectively, the last of which is considered in more 
detail. The modelled time series was constructed without refer- 
ence to any of the composites, except prior to 1990. This means 
that the entirety of cycle 23 is independent of all the TSI com- 
posites. This provides an opportunity to test the model's assump- 
tion: that changes in irradiance are only caused by the evolution 
of surface magnetic flux. 

The SATIRE-S reconstruction is found to agree best with 
the PMOD composite on time scales longer than rotational. On 
shorter time scales the IRMB composite displays a better agree- 
ment, though all three composites are similar on this shorter time 
scale. The model is not just extremely successful at recreating 
rotational and cyclical variability in PMOD, but also in repro- 
ducing the secular trend, particularly apparent in the decline in 
TSI of ~0.2 WirT 2 between the minima of cycle 23. The differ- 
ence between the model and PMOD minima of cycle 21/22 and 
cycle 22/23, relative to the 2008 minima, is less than or equal to 
0.06 WirT 2 and within the uncertainty of the model and of the 
observations. The fact that the model and the PMOD composite 
agree very well provides support for both of them. 

The model is able to account for 92% of the variability in 
PMOD since 1978 and 96% of the variability over cycle 23. 
Due to the general independence of the model, these results sup- 
port the PMOD composite as being the most realistic record of 
TSI. It also suggests that rotational, cyclical and secular vari- 
ations in TSI are mainly the result of the evolution of photo- 
spheric magnetic fields. This finding is in opposition to Frohlich 
(2009| |20n"] > and |Steinhil"b"er[ ( |20 1 0[ ), who proposed that secular 
changes in TSI are caused by a global temperature change (see 
also [Tipping et al. (2007) and references therein). 

It is clear from Fig. 10 that there is disagreement between the 



model and PMOD on periods of one to two years, most signifi- 
cantly over the maximum years of cycle 21 and 22 during which 
KP/512 are employed in the model. KP/512 data are subject to 
many more problems than KP/SPM and SoHO/MDI. These have 
been discussed in sections [3 . 1 ,3| |4.2| and |4.4| 

Some of the differences are also a result of inadequate degra- 
dation corrections in PMOD. During the period prior to 1982 the 
NIMBUS-7/HF radiometer was the only instrument in operation. 
The degradation correction applied to data from this instrument 
was based on the work done on the VIRGO/PM06V radiometer 
in conjunction with ACRIM I data ( |Frohlich||2003| ). While we 
accept that this is a reasonable correction, we believe that some 
uncertainty in PMOD remains during this early period. 

|Krivova et al.| ( |2011a I found that the degradation correction 
of VIRGO during the initial years of observation over 1996-1998 
may not be entirely adequate. Our study adds to that comparison, 
which also used UARS/ACRIM II and DIARAD and PM06V 
radiometers within the SoHO/VIRGO instrument. We find that 
SATIRE-S reconstructions using either SoHO/MDI or KP/SPM 
images disagree in an almost identical manner to the trend in 
PMOD during the first 18 months of the rising phase of cycle 23. 
We did find, however, that the magnetic flux in SoHO/MDI mag- 
netograms is different before and after the period when SoHO 



was not in operation in 1998. This can be attributed to instru- 
ment focus changes. While they can be used in shorter recon- 
struction stretches which agree with the reconstruction using 
KP/SPM (see Fig. BJ, MDI magnetograms should not be used 
to test secular variation, without adequately correcting for this 
effect. 

While TSI is an important observable of changes in solar ir- 
radiance, many recent climate modelling attempts attribute so- 
lar spectral irradiance as physically important when trying to 
understand the impact the Sun has on the Earth's climate. By 
design, SATIRE-S automatically produces spectra over climate- 
relevant wavelengths. A following publication will analyse and 
make these data available. 

The SATIRE-Sind full reconstruction covers 6379 days 
between 1974 and 2009. The daily and smoothed data, 
along with the uncertainty range can be accessed at 
http://www.mps.mpg.de/projects/sun-climate/data.html 
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